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LUNAR-SOLAR PERTURBATIONS 
OF A 24-HOUR SATELLITE 

by 

B. K. Martynenko 


The lunar- solar gravitational perturbations for two hypothetical artificial 
Earth satellites with a sidereal period of revolution of 23^56 m 04 S are analyzed. 


ABSTRACT 

Study of the lunar-solar gravitational perturbations of two hypothetical 

earth satellites having a sidereal period of 23^56 m 04 S . The orbit of the first 
satellite lies in the equatorial plane while that of the second is perpendicular 

h 

to the ecliptic plane. At the initial moment (12 E.T. , January 0, 1950) , the 
nodes of the orbits of the two satellites and the moon coincide. The perigee 
distance of the first satellite is chosen so as to minimize the secular excitations 
of the orbital eccentricity at fixed values of all the remaining elements. The 
orbital parameters of the second satellite are chosen according to the con- 
dition of maximum secular perturbations of the orbital eccentricity. The orbital 
eccentricity of the first satellite is 0. 1, while that of the second is 0. 6. Results 
show that the orbit of the first satellite is stable and orbital eccentricity 
decreases at the rate of 0. 00127 per year. Rapid increase in the eccentricity 
of the second satellite limits its lifetime to 5. 5 years. Secular and periodic 
first order perturbations are calculated. 


INTRODUCTION 

A study is made of the lunar-solar gravitational perturbations of two 
hypothetical artificial Earth satellites having a sidereal period of revolution of 

23^56 m 04 S , i.e., one sidereal day. 

The orbit of one of these satellites (we will call it for reference AES-1) 
lies in the equatorial plane while that of the second (AES-2) lies in a plane 
perpendicular to the ecliptic plane. 
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The eccentricity of the orbit of the equatorial satellite is equal to 0 . 1 
and that of the orbit of the second satellite is equal to 0. 6, i. e. , slightly less 
than Laplace limit. 

h. 

At the initial instant (12 E.T. , January 0, 1950) the nodes of the orbits 
of both satellites and of the Moon coincide while the satellites themselves are 
in the ascending nodes of their orbits. The longitude of the perigee of the AES-1 
was chosen so as to minimize the secular perturbations of the eccentricity of 
its orbit at fixed values of all of the remaining elements. 

The orbital parameters of the AES-2 were chosen in accordance with 
the condition of maximum secular perturbations of the eccentricity of its orbit. 
Being connected with the perigee distance by the relation 

P = a(l - e), 

this element is the most important one in the evolution of the orbit since the 
lifetime of the satellite depends on the character of its variation. 

As shown by the calculations carried out, the orbit of the equatorial 
satellite is stable. Because of a decrease in the eccentricity of its orbit bv 
0. 00127 per year attributable to the combined action of the Moon and the Sun, 
a rise of the perigee altitude takes place at the rate of 53 km per year. In the 
course of time, as the satellite's orbit approaches a circular orbit, the rate 
of this rise decreases. 

A rapid increase in the eccentricity of the orbit of the AES- 2 (by 0. 046 
per year) leads to the termination of the existence of this satellite in only 5. 5 
years. 


The calculation of the lunar-solar perturbations of the first order 
relative to the parameter 


Mo 


= 

m 0 


(m 0 — the Earth's mass; m 2 — the mass of the perturbing body; a — the ratio 
of the major half-axes of the orbits of the satellite and the perturbing body) was 
carried out on the basis of the procedure described in the author's article 
(Martynenko, 1966) on the assumption that the elements of the Moon's orbit 
are constant. 
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The influence of the effects of the noncentricity of the Earth's 
gravitational field, atmospheric resistance, and solar radiation pressure on 
the magnitude of the perturbations attributable to the Moon and the Sun was not 
taken into account. The mutual influence of the lunar and solar perturbations 
also was not taken into consideration. 

The elements of the orbits of the satellites, the Moon, and the Sun 
are shown in Table 1. 


TABLE 1 


THE ELEMENTS OF THE ORBITS OF THE 
AES-1, AES-2, THE MOON AND THE SUN* 


AES-1 

AES-2 

The Moon 

The Sun 

a 

6. 61066874a 

6. 661066874a 

60. 323a 

23370. 1532a 


e 

e 

e 

e 

e 

0. 1 

0. 6 

0. 054900489 

0. 016751 

i 

23f 445788 

o 

o 

o 

5. 0 14538864 

0 

12 

12. 139268 

12. 139268 

12. 139268 

0 

7T 

57. 139268 

147. 139267 

208. 788281 

282. 08039 

€ 

12. 139268 

12. 139268 

’ 57. 786764 

279. 588317 

n 

360. 985649 

360.985649 

13. 1760685 

0.985647344 


* Here a^is the equatorial radius of the Earth. 


1. The First-Order Secular Perturbations 

On the basis of the differential equations for the osculating elements 
of a satellite's orbit (Lyakh, 1960) : 
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x “ m ~ 2 x e s x e ' s x x 

+ s ! 
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q - Ij I = m (mod 2), Ik I = s' (mod 2), 111 = s (mod 2) , 

L = nt + s - t, L' = n't + V - t\ 

t = 12 + N, t' = fi' + N', 

II = 7 T - T, II = 7 T» - T', 

. 1 

sin i sin (12 - 12 ') 


sin N = 
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sin N' 


sin i sin (fl - £2') 
2\/rp 


cos N 


sin i cos i' - cos i sin i' cos - fl') 
2 sfuji 


>T . -sin i' cos i + cos i' sin i cos (fl - fi') 
cos N' = 1 


. 2 I . 2 S2 - O' . ? i + i' 9 S2 — S2 * oi-i’ 

v = sin' - = sin" — sin 2 — - — + cos 2 — sin 2 -^— , 


2 I . ? O - S2' ? i + i' 9 S2 - S2* 9i-i' 

A* = cos"- = sin" — cos" — — + cos 2 — cos 2 , 


a = — , <p = arc sin e, 
a* 


a, e, i, 0, it, and ( are the elements of the satellite's orbit. 


a, e', O', 7 r' , and < ' are elements of the orbit of the perturbing body, 

n and n' are the mean motions of the satellite and the perturbing 
body 

and also on the basis of the tables of the coefficients A m ’. S ’ S , (v , u) and 
(m) q»j»k,l 

a k , n) (Martynenko, 1965) it is easy to obtain expressions for the first- 

order secular perturbations of the elements of the satellite's orbit relative to 
the parameter Indeed, taking into account that the secular terms in the 
perturbations of the elements are generated only by those terms of the 
expansion of the function of R, which explicitly do not contain the time, we obtain 
the following with accuracy to the sixth-order terms relative to the parameters 
a, e, e' 

[6a] = 0; 

[ 6e ] = t x q 2 Yj a m_2 e S e ,S lA m, . S ’ S (r,p) x sin(-kll' - III) 
q +- k= 0 C| ’ •*’ ’ 

j + 1 = 0 

= t X q 2 {-7. 5e 2 ( 1 + 1.5e' 2 + 1. 875 e'Vpsin 211 
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- 3. 75a; ee' (1 + 2. 5e' 2 ) [ (1. 25 m 2 + 2. 5v 2 - 1) M sin (IP - II) 

- (1.25v 2 + 2.5 M 2 “ 1)^ sin (IP + II)] 


- 2. 8125Q!e 3 e'[8. 75 vm 2 sin (IP - 311) 

+ (1.25m 2 + 2.5p 2 - 1) M sin (IP - II) 

- 8.75v 2 Msin (IP + 311) 

- (1.25y 2 + 2.5// - l)v sin (IP 
- 19. 68750! 2 e 2 [ 1. 75 (v 2 + m 2 ) - 1] vm sin 211 
14. 765625Q! 2 e 2 e ,2 


II)] 


(1.75(m 2 + v 2 ) - 1) vm sin 211 

O 


- ( ^ m 2 + 3. 5^ 2 - iJm 2 sin (2IP - 211) 


+ v 2 + 3.5m 2 - ijv 2 sin (2IP + 211) 

57. 6796875 ojWm 2 sin 411 

- 13. 125 o; 3 ee' [(2.625// + 15. 75i^ 2 m 2 

+ 7. 875v 4 - 3.5 m 2 ~ 7y 2 + 1) m sin (IP - II) 

- (2.625i/ + 15. 75y 2 M 2 + 7.875// - 3.5v 2 
-7 // + l)v sin (IP + II) ] 

- 118. 125a/e 2 [4. 125(v 4 + M 4 ) + 12.375 i^ 2 m 2 

- 4. 5 ( v 2 + M 2 ) + 1 ] i^M sin 211 } ; 


t x 03 J a m 2 e S e ' S ' [A m ’. S ’ S ' (p, /u) (lp 2 + kp 3 )sin(-kIP - III) 
q + k-0 q ’ J ’ ’ 

j + 1= 0 


- p 4 ( Av - Am) cos(-kIP - III)] 


= t X qg {-1.5p 4 (l + 1. 5e' 2 + 1. 875e' 4 - 2.1875e’ 6 + 1.5e 2 + 2.25e z e' 2 
+ 2. 8125e 2 e' 4 ) (v - fi) 

- 3. 75e 2 ( 1 + 1. 5e' 2 + 1. 875e' 4 ) [ 2p 2 vju sin 211 + p 4 (M - v) cos 211] 

- 3.75aee'(l + 2. 5e' 2 + 0. 75e 2 ) [(p 2 - p 3 )p ( 1. 25P 2 

+ 2. 5u 2 - 1) sin (II* - II) 

- p 4 ( 5 mm - 3.75 M 2 - 2.5m 2 + 1) cos (IP - II) 

- (p 2 + p 3 )m(1.25v 2 + 2.5jl? - 1) sin (II' + II) 

- p 4 ( 3. 75 m 2 + 2.5m 2 - 5 Pfi - 1) cos (II' + II)] 

- 8. 203125a; e 3 e' [ (3p 2 - p 3 )mm 2 sin (II' - 311) - p 4 M(M - 2v) cos (II’ - 311) 

- ( 3p 2 + p 3 ) m 2 m sin ( 311' + II) + p 4 m(m - 2m) cos (3IP + II)] 

- 3.75a 2 (l + 5e' 2 + 13. 125e' 4 + 5e 2 + 25e 2 e' 2 + 1.875e 4 )p 4 

x ( 1. 75m 3 + 3. 5mm 2 - 3. 5m 2 m - 1.75/^ + n - v) 

- 2. 8125 a 2 e' 2 ( 1 + 3.5e' 2 + 5e 2 ) [ 2pgMM (1. 75m 2 + 1.75m 2 - 1) sin 211' 

+ p 4 ( 5. 25 m 2 m + 1. 75/^ - fi - 1. 75m 3 - 5.25 mm 2 + v) cos 211'] 

- 19. 6875a 2 e 2 (l + 5e' 2 + 0. 5e 2 ) [2p 2 v/n(l. 75v 2 + 1.75 m 2 - 1) sin 211 

+ p 4 (1.75jc? - 5.25 m// 2 + 5.25 m 2 m - M - 1.75m 3 + m) cos 211] 

+ 3. 9375a 2 e 2 e' 2 [2(p 2 - p 3 ) m 2 (2.1875m 2 + 6.5625m 2 - 1.875) sin ( 211' - 211) 

- p 4 p( 13. 125mm - 8.75m 2 - 13.125m 2 + 3.75) cos (2IF - 211) 

- 2(p 2 + p 3 ) v 2 (2. 1875m 2 + 6.5625 M 2 - 1.875) sin (2IP + 211) 

- p 4 v( 8. 75m 2 + 13.125m 2 - 13.125mm - 3.75) cos (2IP + 211) ] 

- 51. 6796875a; 2 e 4 MM[2p 2 MM sin 411 + p 4 (M - v) cos 411] 
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- 13.125o: 3 ee , [(P2 - p 3 )m(2.625 m 4 + 15.75vV + 7.875v 4 

- 3.5 m 2 - 7v 2 + 1) sin (IF - II) 

- p 4 (31.5vM 3 + 31.5 v 3 m - 13.125m 4 - 47.25vV 

- 7.875v 4 - 14i4i 

+ 10.5m 2 + 7v 2 - 1) cos (IF - II) 

- (P 2 + P 3 ) v ( 2. 625v 4 + 15.75vV + 7.875/^ 

- 3.5v 2 -7/^ + 1) sin (IF + II) 

- p 4 ( 13. 125v 4 + 47. 25v 2 M 2 + 7.875 m 4 - 10.5 v 2 - 7/U 2 - 31.5 v 3 m 

- 31.5vj u 3 + 14 vm + 1 ) cos (IF + II)] 

- 6 . 5625o 4 ( 1 + 10. 5e ’ 2 + 10 . 5e 2 )p 4 (4. 125v 5 + 24. 75 v 3 m 2 

+ 12.375vm 4 - 4. 5v 3 

- 9vp 2 + v - 12.375 v 4 m - 24.75 vV - 4.125 M 5 
+ 9v 2 m + 4. 5$ - m) 

- 16. 40625O! 4 e ,2 [ 2p$vn ( 4. 125v 4 + 12.375 v 2 m 2 + 4.125m 4 - 4. 5v 2 

- 4. 5 y^ + 1) sin 2IF 

+ p 4 ( 20 . 625v 4 M + 37.125vV + 4.125^ - 13. 5v 2 m - 4. 5j£ + M 

- 4.125v 5 - 37.125vV - 20.625 vm 4 + 4. 5v 3 
+ 13.5vm 2 - v) cos 2IF] 

-59. 0625 0 4 e 2 [2p 2 pv( 4.125v 4 + 12.375vV + 4.125 m 4 - 4. 5v 2 

- 4. 5m 2 + 1) sin 211 

+ p 4 ( 20. 625v 4 p + 37.125 vV + 4.125 m 5 - 13.5v 2 M 
- 4.5 a? + M" 4.125v 5 
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- 37.125v 3 p 2 - 20.625vm 4 + 4. 5v 3 + 13. 5vp 2 - v) cos 211] 

- 9. 84375q! 6 P4(11. 171875v 7 + 134. 0625vV + 201. 09375vV 

+ 44.6875vp 6 

- 17.875v 5 - 107. 25vV - 53.625 vp 4 + 8. 25v 3 + 16.5vp 2 - v 

- 44.6875vV - 201.09375vV - 134. 0625v 2 jU 5 - 11.171875m 7 
+ 53.625vV + 107.25vV + 17. 875p 5 - 16. 5v 2 m 

- 8. 25ju 3 + ju) }; 


t x q 4 Yj a m " 2 e S e' S [a” 1 ’. 83 ^ < y >M) dPs + sin(-kll' - III) 
q + k= 0 q,J ’ ’ 

3 + 1 = 0 


+ p 7 (Av - A p) cos (-kir - III)] 

t x q 4 { 1 . 5(1 + 1.5e' 2 + 1.875e' 4 - 2.1875e' 6 + 1.5e 2 + 2.25e 2 e' 2 
+ 2. 8125e 2 e' 4 ) x p 7 (v - ju) 

- 3. 75e 2 ( 1 + 1. 5e' 2 + 1. 875e' 4 ) [ 2p 5 vp sin 211 - p 7 (p - v) cos 211] 

-3. 75a;ee'(l + 2.5e' 2 + 0.75e 2 ) [(p 5 - Pg)/4 1.25 m 2 + 2. 5v 2 -l)sin(II' - I 

+ p 7 (5vp - 3.75 m 2 - 2. 5v 2 + 1) cos (II' - II) 

- (P 5 + Pg) v ( 1 . 25v 2 + 2. 5f? - 1) sin (II' + II) 
+ p 7 ( 3.75v 2 + 2.5 m 2 - 5v/u - 1) cos (II' + II)] 
8. 203125 ae 3 e’[ (3p 5 - p 6 )^M 2 sin (II' - 311) + P 7 m(M - 2v) cos (II' - 311 
- (3p^ + p 6 ) v 2 p sin ( 311' + II) - P 7 v(v - 2\x ) cos (311* + II) ] 

+• 3.75q! 2 (1 + 5e' 2 + 13. 125e' 4 + 5e 2 + 25e 2 e' 2 + 1.875e 4 )p 7 
x (1.75v 3 + 3.5 v/U 2 - 3.5v 2 m ~ 1.75^? + n - v) 


- 2. 8125a 2 e' 2 ( 1 + 3.5e' 2 + 5e 2 ) [ 2p 6 ^ ( 1. 75i' 2 + 1.75 a? - D sin 2IP 

- p 7 (5. 25i > 2 n + 1. 75 A? - /U - 1. 75m 3 - 5. 25v/u 2 + v) cos 211* ] 

- 19. 6875o; 2 e 2 ( 1 + 5e f2 + 0. 5e 2 ) [2p 5 t»/u(l. 75i^ 2 + 1.75 a? - 1) sin 211 

- p 7 (1.75a? - 5.25 vfj. 2 + 5.25v 2 /u ~ M ~ 1 . 75m 3 + v) cos 211] 

+ 3. 9375c>! 2 e 2 e' 2 [ 2(p 5 - p G ) a?( 2. 1875/? + 6.5625i, 2 - 1.875) sin(2II r - 211) 
+ p 7 Ai( 13.125i m - 8.75a? - 13.125y 2 + 3.75) cos (2IP - 211) 

- 2(p 5 + p e ) v 2 {2. 1875u 2 + 6.5625a? - 1.875) sin ( 2IP + 211) ] 

+ p 7 y ( 8. 75 p 2 + 13.125a? - 13.125^ - 3.75) cos (2IP + 211) 

- 51. 6796875 o ,2 e 4 i-’Ai 1 2p 5 ^A< sin 411 - p 7 (ju - v) cos 4IIJ 

- 13. 125a?ee' [ ( p 5 - p 6 ) /.t( 2. 625a/ 4 + 15.75i>V + 7.875i^ 4 - 3.5 a? 

- 7j/ 2 + 1) sin (II' - II) 

+ p 7 (31 . 5i'/u 3 + 31 . 5 m 3 /u - 13.125a? - 47. 25m 2 a? 

- 7. 875m 4 - 14un 

+ 10.5a? + 7i^ 2 - 1) cos ( II* - II) 

- (P 5 + p 6 )i'(2.625i> 4 + 15.75i'V + 7.875 a? - 3. 5i' 2 

- 7 At 2 + 1) sin (IP + II) 

+ p 7 (13.125i/ + 47. 25 m 2 a? + 7.875a? - 10.5i> 2 
- 7 A? - 31. 5 m 3 At - 31.5P/I 3 
+ 1 4 m/j -i- 1) cos (IP + II) ] 

+ 6. 5625o?( 1 + 10. 5e’ 2 +10. 5e 2 ) p 7 ( 4. 125i' 5 + 24.75 i, 3 a? 

+ 12. 375^ 4 - 4. 5p* 

- 9ufi 2 + v - 12.375i'V - 24. 75 p 2 a? - 4.125a? 
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+ 9 m 2 m + 4. 5a? - A ( ) 


16. 40625a' 4 e' 2 [2p 6 MM(4. 125 m 4 + 12.375mV + 4.125/u 4 - 4.5 m 2 

- 4. 5 jU 2 + 1) sin 2IP 

- p 7 (20.625fV + 37.125m 2 a? + 4.125a? - 13. 5m 2 p - 4. 5a? + p 

- 4.125m 5 - 37.125 mV - 20. 625m/u 4 + 4.5m 3 

+ 13.5M/U 2 - m) cos2II’] 

59. 0625a 4 e 2 [2p 5 /iM(4. 125 m 4 + 12. 375m 2 A? + 4. 125p 4 - 4.5 m 2 

- 4.5a? + 1) sin 211 

- p 7 (20. 625 m 4 ju + 37.125mV + 4.125a? - 13.5m 2 p 

- 4. 5a? + p 

- 4.125m 5 - 37.125 mV - 20. 625m/u 4 + 4.5m 3 

+ 13.5mp 2 - m) cos 211] 

+ 9. 84375o' 6 p 7 ( 11. 171875 m 7 + 134. 0625m 5 A? + 201. 09375m V 
+ 44. 6875mm G 

- 17.875m 5 - 107. 25m 3 a? - 53.625mp 4 + 8.25m 3 + 16.5 m/u 2 - m 

- 44. 6875m g p - 201. 09375mV - 134. 0625m 2 a? - 11.171875a? 

: 53. 625m 4 p - 107. 25 m 2 a? + 17. 875a? - 16.5m 2 p 

- 8. 25a? + P) }; 

. , y m - 2 s - 1 s' 

q 5 [ ] -r t • q 2 h scv e e' 

q * k = 0 

j : 1-0 

. m, s, s' , . , ... 

A . . , (m , p) cos ( — kll* - ill) 

q,J,K. 1 
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= q 5 [ SS2 ] + t x q 2 {l. 5e[5(l + 1. 5e' 2 + 1. 875e’ 4 ) cm cos 211 
+ (1 + 1. 5e’ 2 + 1. 875e' 4 ) (1. 5c 2 + 1.5/? - 1) ] 

- 3.75ae'((l + 2.5e' 2 + 2. 25e 2 ) ( 1. 25// + 2.5c 2 
- 1) m cos (IP - II) 

+ ( 1 + 2. 5e’ 2 + 2. 25e z ) ( 1 . 25c 2 + 2.5// - 1) c cos ( IP + II) ] 
- 24.609375ae 2 e’ [ m cos (IP - 311) + c cos (IP + 311) ] cm 
+ 3. 75o/e[ ( 1 + 5e' 2 + 0. 75e 2 ) ( 4. 375c 4 + 17. 5 v 2 ju 2 
+ 4. 375 m 4 - 5c 2 - 5// + 1) 

+ 10.5(1 + 5e' 2 + e 2 ) (1.75c 2 +■ 1.75m 2 - 1 ) cm cos 211] 


+ o/ee' 2 


14.765625 ( ^ // + 3.5c 2 - l)// cos ( 2IP - 211) 

\ J V 


-) 


+ 28. 125(1. 75c 2 + 1.75// - 1) vfi cos 2IP 


( 


105 2 o c 2 

— c 2 + 3.5 m 2 


- l^c 2 


cos ( 2IP + 211) 


+ 14.765625 
+ 103. 359375 o/eV// cos 411 

- 13. 125c/e'[ (2. 625/u 4 t- 15. 75c 2 ./ ‘ 7.875c 4 - 3.5// 
- 7c 2 + 1) n cos ( IP - II) 


+ ( 2. 625m 4 + 15. 75c 2 m 2 + 7.875/u 4 - 3.5c 2 
- 7// 1) c cos (IP t- II) ) 

+ 6. 5625a 4 e| 18(4. 125m 4 + 12. 375c 2 // f 4.125c 4 - 4.5// 

- 4. 5// f- 1) i-’M cos 211 

14.4375c 6 ■ 129.9375cV ' 129.9375c 2 m 4 14. 4375m® 

- 23.625c 4 - 94. 5cV - 23.625m 4 f 10.5c 2 : 10.5m 2 - 1 J > ; 
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[6c] = cfc!^] + t x 


m — 2 s - 1 .s' 

a ( sp 0 - 2me) e e' 


= q 5 [« J ] 


- :s. 


+ 3. 


4. Z 

q + k= 0 
j ^ 1=0 

A m,S l S \ (v,n) cos (-kll* - HI) 

q, j,k, 1 

t t x q G { [ 1. 5e(p 0 - 2e)(l + 1.5e' 2 + 1 . 875e' ) 2 - 3e 

3. 75e' 4 - 4.375e' 6 ] 

x ( 1 . 5v 2 + 1.5 1 ? - 1) + 7.5e(p 0 - 2e) (1 + 1.5e' 2 
i 1. 875e' 4 ) via cos 211 

75ae'[(<Pt> - 6e)(l + 2.5e' 2 ) + 2.25 <p 0 - 2e) e 2 ) ( ( 1. 25p 2 + 2. 5^ 2 
- 1) IX cos ( II’ - II) 

+ (1.25V* t- 2.5 m 2 - 1)m cos (II’ + II)) 

+ 6. 5625 ( p„ - 2 e)eV (mcos (IP - 311) + v cos (II' + 311))] 
75o? 2 |(-- 0.8 - 4e’ 2 - 10. 5e' 4 + e(p 0 - 4e)(l +- 5e' 2 ) 

+ 0. 75c 3 (p () - 2e)) 

x ( 4. 375v 4 i- 17. 5v 2 p 2 ► 4.375 m 4 - 5^ 2 - 5p 2 i- 1) 

+ 10. 5e ( ( p„ - 4c) (1 + 5e' 2 ) + (p„ - 2e) e 2 ) ( 1. 75 m 2 
+ 1. 75i> 2 - 1) I'M cos 211 

h e »2 (_ 4 . 5 _ 21e' 2 + 7.5e(p 0 - 4e)) (1.75v 2 + 1.75m 2 
- 1) VH COS 211' ] 

. lee’ 2 (p 0 - 4e) | (2. 1875M 2 + 6.5625i' 2 - 1. 875) m 2 cos (2IP - 211) 

+- ( 2. 1875y 2 + 6.5625 m 2 - 1. 875) v* cos (2IP + 211) ] 

+ 27. 5625e 3 (p 0 - 2e) v 2 tf cos 411] 
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- 13. 125o; 3 e , (Po - 10e) [ (2. 625ju 4 + 15. 75vV + 7.875v 4 - 3.5m 2 


- 7 v 2 + 1) m cos (II» - II) 

+ (2.625v 4 + 15. 75 v 2 m 2 + 7. 875M 4 - 3. 5v 2 

- 7m 2 + 1) v cos (IP + II) ] 

- 3.75Q! 4 [(l + 10. 5e' 2 - 1.75e(p 0 - 6e)) (14.4375V 6 + 129.9375vV 

+ 129. 9375v 2 M 4 

+ 14.4375/u 6 - 23.625V 4 - 94.5v 2 M 2 - 23.625M 4 
+ 10.5v 2 + 10. 5 M 2 - 1) 

+ (52. 5e’ 2 cos 2IP - 31.5e(p 0 - 6e) cos 2Il) (4. 125v 4 
+ 12. 375 v 2 m 2 

+ 4. 125M 4 - 4. 5v 2 - 4. 5 m 2 + l)vju] 

- 16 a 6 ( 13. 746643v 8 + 219. 946289v 6 M 2 + 494. 87915vV 

+ 219. 946289v 2 m 6 

+ 13.746643m 8 - 29.3261718v 6 - 263. 935546v 4 M 2 

- 263. 935546v 2 

- 29.3261718m 6 + 20. 3027343v 4 + 81. 2109375 v 2 m 2 

+ 20.3027343m 4 

- 4.92l875v 2 - 4.921875 m 2 + 0.2734375)}; 

[fin] = 0. 

Here 

q 2 = m o n - e 2 /e; p 0 = e n/i - e 2 /(l + Nl - e 2 ) ; 

qa = t^n/41 " e 2 ; P 2 = cos N (P - v)/2n/vm"; 
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Q4 = 

A^n/ (n/i - e 2 sin i) ; 

P 3 = 

cos N/Z’Jw; 

Q5 = 

tan ^ sin i; 

P4 = 

^ sin N; 

Qe = 


P5 = 

sin N ( p - v) /2N r i 



Pe = 

sin N/ 2 N/i 4 t; 



P? = 

N/vju cos N. 


The values of the parameters contained in these expressions are shown 
in Table 2. 


TABLE 2 


ELEMENTS OF THE RELATIVE POSITION OF THE 
ORBITS OF THE SATELLITES, THE MOON AND THE SUN 


Per- 

AES-1 

AES-2 

turbing 

The 

The 

The 

The 

Body 

Moon 

Sun 

Moon 

Sun 

a 

0. 109587864 

0. 000282868 

0.109587864 

0. 000282868 


0.000016173 

0. 000007455 

0.000016173 

0.000007455 

V 

0. 025288355 

0. 041281524 

0. 455158342 

0.5 

V 

0.974711645 

0.958718475 

0.544841658 

0.5 

N 

0 

0 

0 

0 

N' 

0 

12.°139268 

0 

12.°139268 

t ,r T 

12.° 139268 

12T139268 

12.°139268 

12.°139268 

ii 

45° 

45° 

135° 

135° 

ip 

196°649013 

269.°941122 

196.°649013 

269.°941122 

Po 

0. 049874371 

0. 049874371 

0.26666666 

0. 26666666 

P2 

3.02364997 

2. 3058063 

0. 090046172 

0 

P3 

3.18472277 

2. 51331309 

1. 00404597 

1 

P4.P5.P6 

0 

0 

0 

0 

P? 

0. 156999536 

0. 198940594 

0. 497985166 

0. 5 

Q2 

0. 0010138 

0. 000467355 

0. 000135863 

0. 000062628 

Q3 

0. 0001024 

0. 000047208 

0. 000127372 

0. 000058714 

Q4 

0. 00025739 

0. 000118647 

0. 000127372 

0. 000058714 

P5 

0. 082563048 

0. 082563048 

1 

1 

q6 

0. 000101897 

0. 000046971 

0. 000101897 

0.000046971 
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Table 3 contains numerical values of the secular changes in the elements 
of the orbits of the satellites AES-1 and AES-2 attributable to the Moon, the 
Sun, and the Earth's oblateness after 1 year. The last were calculated by the 
following formulas (Proskurin, Batrakov, 1960): 


[S2] = -J 


a 


a 


cos 1 


(1 - e 2 ) 


2. 2 ^ ^ 365. 25, 


[7r] = 2 J 


5 cos 2 i - 2 cos i - 1 


(1 - e 2 ) 


2 \ 2 


n x 365. 25, 


r , t , j e 1 (5 + 3 \l 1 - e 2 ) cos 2 i - 2 cos i - 1 - 

|£i “ 2 J \T/ (1 - e’) ! 


n x 365.25 


with a e = 6, 378. 39 km and J = 0. 00164147. 


TABLE 3 


SECULAR PERTURBATIONS OF THE ELEMENTS OF THE 
AES-1 AND AES-2 ORBITS ATTRIBUTABLE TO THE 
MOON, THE SUN, AND THE EARTH'S OBLATENESS 



[a] 

[e] 

[i] 

[Q] 

[7T] 

[e] 

( The Moon 

0 

-0. 00075 

-0.°016 

-If 3 

l.°3 

-If 9 

AES-1 1 The Sun 

0 

-0. 00052 

-0. 007 

-0. 69 

0. 5 

-0. 8 

( The Oblateness 

0 

0 

0 

-5 

4.7 

10 

( The Moon 

0 

0. 030 

0. 15 

-0. 28 

-0.9 

1 

AES-2 } The Sun 

0 

0.016 

0 

0 

-0.3 

0. 7 

v The Oblateness 

0 

0 

0 

4.7 

3. 6 

1. 1 


In these formulas "i" is the angle of inclination of the satellite's orbit 
to the plane of the equator. 

As may be seen from Table 3, on the average the secular perturbations 
of the elements of the satellites' orbits attributable to the Moon are two to three 
times greater than those attributable to the Sun. 

In their turn, the aggregate values of the lunar-solar perturbations 
of the elements are two to three times smaller than the respective values 
explained by the Earth's oblateness. 
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In comparing the secular lunar-solar perturbations of the AES-1 with 
those of the AES-2, in addition to the difference in the magnitude of the 
perturbations of the eccentricities already mentioned in the Introduction, 
attention is drawn by the considerable preponderance of the secular perturbations 
of the inclination of the AES-2 orbit (because of the Moon) over the respective 
magnitude for the AES-1. The secular changes in the elements fi, ir, and e 
caused by the influence of the Moon and the Sun are several times smaller for 
the AES-2 than for the AES-1. 

It is interesting to note that the elements "i" and S2 of the AES-2 orbit 
undergo practically no secular perturbations from the Sun. This is explained 
by the fact that P2 = P4 - P5 = Pe = v - n = 0. Because of this, the greater 
portion of the addends in the expressions for the secular perturbations of these 
elements vanish. 


2. The First-Order Periodic Perturbations 

Graphs shown in Figures 1-24 (Appendix) were constructed using the 
points with the time intervals of 1 day for the lunar perturbations and 10 days 
for the solar perturbations in order to secure a general idea of the character 
of the changes in the elements of the orbits of the satellites being examined, 
under the effect of the gravitational perturbations attributable to the Moon and 
the Sun. In addition to evaluating the magnitude of the secular perturbations 
it was intended to secure information concerning the periodic changes in the 
elements whose period is comparable in length with the duration of one 
revolution of the perturbing body. Therefore, the graphs constructed cover 
the time intervals of 70 days for the lunar perturbations and 700 days for the 
solar perturbations . 

In Figure 25 ( Appendix) are shown the secular changes in the elements 
of the AES-2 orbit during 17 years. 

We will examine some of the special characteristics of the periodic 
perturbations of the elements of the AES-1 and AES-2 attributable to the Moon 
and the Sun. 


AES-1 

Perturbations of the major half-axis: 


6 a 


max 


Lunar — the period is about 30.2 days; 6a . = -1.2 km, 

min 

= 2.3 km; 
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Solar — the period is 1 year; during this time, two minima of different 
level are observed: -90 m and -360 m, and two maxima: -40 m and -45m; the 
minima correspond to the angular departures of the Sun from the perigee by 
90 and 270°. 

Perturbations of the eccentricity: 

Lunar — the period is equal to the sidereal lunar month (27-1/3 days) ; 
the amplitude is about 20 x 1 0 — 5 ; 

Solar — the period is 1 year; the amplitude is 42 x 10 -5 ; two minima 
and two maxima during the period set in respectively when the Sun is at 45, 225, 
135, and 315° from the perigee. 

Perturbations of the inclination: 

Lunar — the period is 27-1/3 days; the amplitude is Of 003; 

Solar — the period is 1 year; the amplitude is Of 023; the minima 
correspond to the angular departures of the Sun from the perigee by 0 and 180°; 
the maxima set in when u = 90 and 270° (u is the true anomaly of the Sun) . 

Perturbations of the longitude of the ascending node: 

Lunar — the period is 27-1/3 days; the amplitude is Of 006; 

Solar — the period is 1 year; the amplitude is Of 05; the 6ft curve has 
two inflection points corresponding to u = 90 and 270°. 

Perturbations of the longitude of the perigee: 

Lunar — the period is 27-1/3 days; the amplitude is Of 13; 

Solar — the period is 1 year; the amplitude is Of 38; the maxima are at 
v = 0 and 180°, the minima at u = 90 and 270°. 

Perturbations of the mean longitude per epoch: 

Lunar — the period is 27-1/3 days; the amplitude is 0.°0035; 

Solar — the period is 1 year; the amplitude is Of 015; the maxima 
are at u = 60 and 240°, the minima at u = 150 and 330°. 
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Thus, clearly marked periodic components whose periods coincide with 
the sidereal periods of the revolution of the perturbing bodies are observed in 
the osculating elements of the AES-1 orbit in the time intervals examined. 

An exception is that of lunar perturbations of the major half-axis whose 
period is equal to 30. 2 days. 

And what is more, the extremal values of the solar perturbations of 
the elements regularly depend on the position of the Sun in relation to the perigee. 

The periodic solar perturbations of the elements of the AES-1 orbit are 
2 to 10 times greater than the perturbations of the respective elements 
attributable to the Moon. An exception is that of the perturbations of the major 
half-axis, which are 6. 5 times greater from the Moon than from the Sun. 


AES-2 

Perturbations of the major half-axis: 

Lunar — the period is 27-1/3 days; 6a . =-1.7 km; 6a =3 1cm; 

min max 

Solar — the period is 1 year; during the period, three minima are 
observed : 


at u 
at v 
at u 

and three maxima: 

at u 
at u 
at u 


0 

6a = 

0 , 

120" 

6a = 

0 , 

240 c 

6a = 

-0. 5 1cm 


60° 

6a = 

1.6, 

160° 

oa = 

0.95, 

320° 

6a = 

0. 4 1cm. 


Perturbations of the eccentricity: 


Lunar — the period is 27-1/3 days; the amplitude is Of 00025; 

Solar -- the period is 1 year; the amplitude is Of 001; the maxima are 
at u = 45 and 225°, the minima at e = 135 and 315°. 
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Perturbations of the inclination: 


Lunar — the period is 27-1/3 days; the amplitude is Of 01 7; 

Solar — the period is 1 year; during this period two equal minima of 
the 6i = 0 are observed at u = 0 and 180° , and two equal maxima of the 
<5i = Of 23 at u = 90 and 270°. 

Perturbations of the longitude of the ascending node: 

Lunar — the period is 27-1/3 days; amplitude is Of 0115; 

Solar — the period is 1 year; two equal maxima of the = 0 correspond 
to a u = 0 and 180°; two equal minima of the <5t2 = 0. 135 at u = 90 and 270°. 

Perturbations of the longitude of the perigee: 

Lunar — the period is 27-1/3 days; the amplitude is Of 015; 

Solar — the period is 1 year; the amplitude is Of 095; the minima are 
at v = 70 and 250°, the maxima at u = 160 and 340°. 

Perturbations of the mean longitude per epoch: 

Lunar — the period is 27-1/3 days; the amplitude is 0. 025; 

Solar — the period is 1 year; the amplitude is Of 17; the maxima are 
at u = 45 and 225°, the minima at v = 120 and 300°. 

Thus, for the AES-2, as for the AES-1, in the time intervals of several 
revolutions of the Moon and the Sun the main role among all periodic perturbations 
is played by the perturbations whose periods coincide with the sidereal periods 
of the revolutions of these bodies. 

The points corresponding to the extremal values of the solar pertur- 
bations of the elements of the AES-2 are regularly located relative to the 
perigee of the Sun's orbit. 

With the exception of the major half- axis all elements of the AES-2 
orbit undergo periodic perturbations attributable to the Sun that are several 
times greater than those attributable to the Moon. 
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Appendix 


THE LUNAR-SOLAR GRAVITATIONAL PERTURBATIONS OF TWO 
HYPOTHETICAL ARTIFICAL EARTH SATELLITES (THE CURVES) 


FIGURE 1. LUNAR PERTURBA- 
TIONS OF THE MAJOR 
HALF-AXIS OF THE AES-2 
ORBIT. 




FIGURE 2. SOLAR PERTURBA- 
TIONS OF THE MAJOR 
HALF-AXIS OF THE AES-1 
ORBITS. 


FIGURE 3. LUNAR PERTURBA- 
TIONS OF THE ECCENTRICITY 
OF THE AES-1 ORBIT. 


6 a - 10 5 


DAYS 




FIGURE 4. SOLAR PERTURBA- 
TIONS OF THE ECCENTRICITY 
OF THE AES-1 ORBIT. 


FIGURE 5. LUNAR PERTURBA- 
TIONS OF THE INCLINATION 
OF THE AES-1 ORBIT. 
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FIGURE 6. SOLAR PERTURBA 
TIONS OF THE INCLINATION 
OF THE AES-1 ORBIT. 


FIGURE 7. LUNAR PERTURBA- 
TIONS OF THE LONGITUDE 
OF THE ASCENDING NODE 
OF THE AES-1 ORBIT. 


FIGURE 8. SOLAR PERTURBA 
TIONS OF THE LONGITUDE 
OF THE ASCENDING NODE 
OF THE AES-1 ORBIT. 


FIGURE 9. LUNAR PERTURBA- 
TIONS OF THE LONGITUDE 
OF THE PERIGEE OF THE 
AES-1 ORBIT. 


FIGURE 10. SOLAR PERTURBA 
TIONS OF THE LONGITUDE 
OF THE PERIGEE OF THE 
AES-1 ORBIT. 


FIGURE 11. LUNAR PERTURBA- 
TIONS OF THE MEAN 
LONGITUDE OF THE AES-1 
PER EPOCH. 






6Q = io 2 
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FIGURE 12. SOLAR PERTURBA 
TIONS OF THE MEAN 
LATITUDE OF THE AES-1 
PER EPOCH. 


FIGURE 13. LUNAR PERTURBA- 
TIONS OF THE MAJOR 
HALF- AXIS OF THE AES-2 
ORBIT. 


FIGURE 14. SOLAR PERTURBA 
TIONS OF THE MAJOR 
HALF-AXIS OF THE AES-2 
ORBIT. 


FIGURE 15. LUNAR PERTURBA- 
TIONS OF THE ECCENTRICITY 
OF THE AES-2 ORBIT. 






FIGURE 16. SOLAR PERTURBA- 
TIONS OF THE ECCENTRICITY 
OF THE AES-2 ORBIT. 




FIGURE 17. LUNAR PERTURBA- 
TIONS OF THE INCLINATION 
OF THE AES-2 ORBIT. 
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6 i = i o 2 



DAYS 


FIGURE 18. SOLAR PERTURBA- 
TIONS OF THE INCLINATION 
OF THE AES-2 ORBIT. 


FIGURE 19. LUNAR PERTURBA- 
TIONS OF THE ASCENDING 
NODE OF THE AES-2 ORBIT. 
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FIGURE 20. SOLAR PERTURBA- 
TIONS OF THE LONGITUDE 
OF THE ASCENDING NODE 
OF THE AES-2 ORBIT. 


FIGURE 21. LUNAR PERTURBA- 
TIONS OF THE LONGITUDE 
OF THE PERIGEE OF THE 
AES-2 ORBIT. 


6 IT = 10 3 



6 ~ io 2 



- 40k 


- 60 |- 



FIGURE 22. SOLAR PERTURBA- 
TIONS OF THE LONGITUDE 
OF THE PERIGEE OF THE 
AES-2 ORBIT. 
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FIGURE 23. LUNAR PERTURBA- 
TIONS OF THE LONGITUDE 
OF THE AES-2 PER EPOCH. 



FIGURE 25. SECULAR SOLAR 
PERTURBATIONS OF THE 
ELEMENTS OF THE AES-2 
ORBIT: ECCENTRICITY, 
THE LONGITUDE OF THE 
PERIGEE AND OF THE 
MEAN LATITUDE PER 
EPOCH. 

A Stt° x 10 IS LAID OFF ON 
THE LOWER STRAIGHT LINE 
( THE DOT-AND-DASH LINE) . 



DAYS 


FIGURE 24. SOLAR PERTURBA- 
TIONS OF THE MEAN 
LONGITUDE OF THE AES-2 
PER EPOCH. 
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increase in the eccentricity of the second satellite limits its lifetime to 5.5 years. Secular 
and periodic first order perturbations are calculated. 
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